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We study inelastic decay of bosonic excitations in a Luttinger liquid. In a model with linear 
excitation spectrum the decay rate diverges. We show that this difficulty is resolved when the inter- 
action between constituent particles is strong, and the excitation spectrum is nonlinear. Although 
at low energies the nonlinearity is weak, it regularizes the divergence in the decay rate. We develop a 
theoretical description of the approach of the system to thermal equilibrium. The typical relaxation 
rate scales as the fifth power of temperature. 
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One-dimensional interacting systems [T] are funda- 
mentally different from their higher-dimensional counter- 
parts P]. Regardless of the statistics of the constituent 
particles, elementary excitations in one dimension are be- 
lieved to be bosons |1I3I4) . the waves of density. Similar 
to sound waves in ordinary fluids, bosonic excitations in 
such a Luttinger liquid [3 4 have linear spectrum at low 
energies ujq = s\q\. Here q is the wave number and s is 
the velocity. 

Just as quasiparticles in the Fermi liquid [2] , bosons in 
the Luttinger liquid do not represent exact eigenstates of 
a generic one-dimensional system. At finite energies, the 
corresponding effective Hamiltonian should be amended 
by irrelevant in the renormalization group sense pertur- 
bations [3], such as interaction between the bosons. How- 
ever, a naive attempt to account for this interaction per- 
turbatively immediately leads to difficulties. 

Consider, for example, the interaction-induced decay 
of a boson with wave number q into two bosons with 
wave numbers q'^ and q2, see Fig.jlja). The corresponding 
inelastic scattering rate is given by the Fermi golden rule, 

T^^(X Jdq[dq!2[...]5{q-q[-q'2)S{ujq-ujq,^~ujg,J, (1) 

where the two (5-functions express the momentum and 
energy conservation. When all three wave numbers have 
the same sign, the second (5-function reduces to s^^S{q~ 
i'l ~ 92); ^iid the rate ([!]) diverges. 

One way around the failure of the perturbation theory 
is to abandon the effective Luttinger liquid description 
altogether and approach the problem from the original 
fermionic perspective |5l6j . Indeed, for noninteracting 
fermions the spectral weight of the dynamic structure 
factor (Fourier transform of the density-density correla- 
tion function) at a fixed q is spread uniformly over a 
narrow interval of the width 

5ujq — hp^q^/m^ (2) 

about uj — LOq. Here is the effective mass, which for 
free fermions coincides with the bare mass m, and q is the 
dimensionless (measured in units of the particle density 



p) wave number. At sufficiently small g, Eq. ^ is appli- 
cable to interacting fermions as well [5H2 • The inverse of 
the width, l/Sujq^ provides a natural estimate of the life- 
time of bosons in the Luttinger liquid. Since 5ujq oc w^, 
the bosons indeed represent well-defined quasiparticles. 

In this Letter we develop an alternative approach, 
based on the observation that divergences that plague the 
evaluation of the quasiparticle decay rate in the conven- 
tional Luttinger liquid theory can be cured if the boson 
spectrum is nonlinear, such as 

uJq^s\q\{l~S,q^). (3) 

Even for a weak nonlinearity <C 1, decay of a single 
boson into two is forbidden by the momentum and en- 
ergy conservation laws and can only occur virtually. The 
simplest real scattering process involves two bosons both 
in the initial and in the final states, see Fig. [ijb), and 
has a finite rate. 

Keeping the nonlinear correction in Eq. ([3]) is justi- 
fied only in the limit of strong repulsion, i.e., when the 
Luttinger liquid parameter [I] if = TrHp^/ms is small. 
Indeed, the correction must exceed the width Sujq [see 
Eq. ([2])], which can be viewed as an uncertainty in the 
energy of the Luttinger liquid's boson. Using the esti- 
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FIG. 1: (a) For bosons with a linear spectrum scattering 
of a single boson (filled circle) into two (open circles) has a 
divergent rate, (b) For bosons with a nonlinear spectrum the 
simplest scattering event satisfying the momentum and en- 
ergy conservation laws involves two bosons both in the initial 
state (filled circles) and in the final state (open circles). For 
given qi and 52 , the conservation laws yield a unique set , ^2 > 
thus leading to a finite transition rate. 
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mate [S] mt^jm ^ vK^ valid for if <C 1, we arrive at the 
condition » ^/K. 

For if <C 1, Eq. ([3| is applicable in a broad range of 
wave numbers K <^ S^q <^ and spectrum nonlin- 
earity has a dramatic effect on inelastic scattering. For 
the scattering process with two bosons [qi and 52) in the 
initial state and two bosons [q'l and q'2) in the final state 
[see Fig.[IJb)], the conservation laws qi+q2 = 9i + 92 s-^d 
yield a unique set q'l , (72 for given 



^91 



"92 



-92 



qi,q2- Moreover, if qi,q'i, and belong to the same (say, 
right-moving) branch of the spectrum [see Fig. m[b)] , the 
remaining wave number is given hy q2 ~ — {3(,/2)qiq[q'2, 
i.e., the sign of 92 is opposite to that of qi,qi,q2 and 
the momentum transferred from the left-moving branch 
of the spectrum in each act of scattering is parametri- 
cally small compared with that redistributed among the 
three right-moving bosons. Accordingly, the process re- 
sembles decay of a single right-moving boson into two. 
However, unlike for bosons with strictly linear spectrum, 
the mere presence of the left-moving boson with a very 
small momentum, as required by the conservation laws, 
is sufficient to regularize the divergences. 

We describe our strongly interacting system by the 
Hamiltonian 



(4) 



where pi and xi are, respectively, the momentum and 
position of the Zth particle {I = 1, . . . ,iV), and V{x) is 
the interaction potential. In the strong repulsion limit 



(i.e., for d^V/dx^ 



--IIP 



3> h p /m, which is equivalent to 



<C 1) the particles, regardless of their statistics, form 
at low energies a periodic chain, the so-called Wigner 
crystal (see [9] for a review). Although in one dimen- 
sion quantum fluctuations destroy the true long-range 
order [TD] , the interparticle distance remains close to 

Similar to ordinary crystals, the elementary excita- 
tions of the Wigner crystal are phonons. These phonons 
are nothing but the waves of density, with a typical for 
phonons linear dispersion at small momenta, i.e., the 
phonons coincide with the bosons of the effective Lut- 
tinger liquid theory. The boson spectrum Wq in the lead- 
ing (zero) order in H can be found by expanding the 
potential energy in Eq. Q to second order in the dis- 
placements of the particles from the corresponding lat- 
tice sites ui = xi — l/p, and solving classical equations of 
motion [TT]. For small q, this yields Eq. (|3| with model- 
dependent s and ^ . 

Interaction between the bosons arises from the higher- 
order (anharmonic) terms in the expansion of the poten- 
tial energy in Eq. ^ in the displacements ui . A scatter- 
ing process with two bosons both in the initial and in the 
final states, see Fig. ijb), can occur either in the first or- 
der in the quartic anharmonicity, or in the second order in 
the cubic anharmonicity , and the corresponding con- 
tributions to the on-shell scattering amplitude iq-^q^;q'^q'^ 
are of the same order of magnitude. If all four wave num- 



bers are small, the amplitude simplifies [T^] to 



9i92;9i92 - N m 



I / ni/2 
9i929i92 



(5) 



This expression is easy to understand if one notices that 
each boson with wave number q participating in scatter- 
ing contributes a factor of {fi/u:qY/'^\q\ cx [^q])^^"^ to the 
amplitude. The dimensionless parameter A in Eq. ([s]) de- 
pends on the functional form of V{x) |12j . In particular, 
A = for V{x) (X l/sinh^(c/9x) and V[x) cx l/x^ [12], as 
expected for integrable models [T3] exhibiting no relax- 
ation. For a generic interaction potential |A| is of order 
unity. In particular, A = —3/4 for screened Coulomb 
interaction, see [T2j . 

Inelastic scattering leads to the relaxation of the bo- 
son distribution function Nq towards equilibrium. The 
evolution of Nq is described by the Boltzmann equation, 
which for a homogeneous system in the absence of exter- 
nal fields has the form [T3] 



dNq 
dt 



Iout[Nq] +I;n[Nq], 



(6) 



where the two terms in the right-hand side describe, re- 
spectively, the scattering out of single-boson state q, and 
the scattering into this state. In the leading order in /i, 
these terms are given by 



lout [Nq] = Wq,p,q,^q,NqN,il + Nq,){l + Nq,) , 

p 9i>g2 

{l + Nq){l + N,)Nq^Nq, 



P 9l>92 



with 



9l.92;9i.92 



^2 I '?l92;9l92l "91+92^91+92 
X S{ujq^ +UJq, ~UJq, 



92 



)■ (7) 



We begin the analysis of Eqs. (|6|-(7| by consider- 
ing the relaxation rate of a single high-energy boson. 
Specifically, we assume that the distribution function Nq 
differs from its equilibrium form, the Bose distribution 
Hq — (^e''"'/"^ — 1) , in the population of a single state 
with q in the range T/hs q In this limit 

lin [Nq] is exponentially suppressed, and Eq. ^ reduces 
to dNq/dt — —Nq/Tq with the relaxation rate 



487r3 



{T/hs)q\ q « {T/hsS,Y'^ 



a{T/hs 



q » 



(8) 



where a = 32C(3)/3. Here Ci^) is the Ricmann's zeta- 
function, C(3) « 1.2. Although Eq. ^ is not directly 
applicable to thermal bosons with energy of the order of 
T, setting q ^ T/hs in Eq. ([s]) results in a correct order- 
of-magnitude estimate of the typical scattering rate, see 
Eq. (fl6|) below. 
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Independently of the initial state, at i — cxd the distri- 
bution function Nq relaxes to Uq. In order to study the 
approach to equilibrium, we substitute 



Nq^nq+ gqfq, Qq = \J nq{l + Uq) 



(9) 



into Eqs. ([6])-([7]), neglect all but linear in fq contribu- 
tions, and obtain 



dt 



27r 



|2 lg_ ^ fp_ _ fii _ fq2 

.9q 9p 9qi 9q2 
X 9p9qi9q2^i^q+^P^^qi ~ ^^92 ) (^g+p,?! +g2 ' (10) 



fl2 ^ |''«P;9l92 I 

V qi>q2 



The linearized Boltzmann equation ( 10 1 is applicable 



for both positive and negative q. Focusing from now on 
on g > 0, we note that Eq. ( 10 ) simplifies considerably if 



(11) 



The first inequality in Eq. ( 11 ) ensures that contributions 
from the processes with all bosons but q on the left- 
moving branch of the spectrum are exponentially sup- 
pressed. The second inequality in Eq. (Ill guarantees 



that the wave number of the only left-moving boson par- 
ticipating in the remaining scattering processes is much 
smaller than T/Hs. Under these conditions, the spectrum 
in the right-hand side of Eq. ( 10 1 can be linearized, which 
amounts to neglecting corrections of order £^{T/tis)^ <C 1 
[this inequality is implicit in Eq. (11)]. This approxima- 
tion corresponds to the substitution into Eq. ( 10 1 



S{ujq +U}p- UJq^ - UJq^) Sq+p^g^+q^ « ^ [S{p + 0) Sq,q^+q.^ + 5{q2 + 0) (5q+p,gJ , (12) 

where 5{k + {)) indicates that k is an infinitesimal wave number on the left- moving branch. This yields 



^ = - -r^>'^K^s{T/hs)q [ dqi\}-[ dq2 5{q - qi - q2) qigqiq29q2\ — - — - — 

Ot 47r-^ Jo \^Jo \9q 9qi 9q2 



+ 1 dpS{q+p-q^)pgpq^gq,[^ + ^^^]}, (13) 

^ 9q 9p 9qi . 



where 9q,9p,9qi, and gq^_ are given by Eq. ^ with a 
linearized spectrum, e.g., gq = [2 sinh(/is(7/2r)] . The 
factor T/hsin the right-hand side of Eq. ( 13 ) is a remnant 
of the left-moving boson. Indeed, its wave number k [k 
is either p or q2, see Eq. ([T2|] appears in Eq. (fTol) in 



combination \k\gk, where the factor |fc| comes from the 
square of the amplitude ([5|. For |fc| <C T/hs, we have 
= {T/hs)\k\-^, which gives \k\gu = T/hs. 

Note that the parameter ^ [see Eq. ([3|] does not appear 
explicitly in Eq. (13 1. This is consistent with the above 



result for the relaxation rate of high-energy bosons: 

is independent of f at g < (T/fof)^, see Eq. Note 



also that all wave numbers in Eq. (13) are strictly posi- 
tive: coupling between bosons moving in opposite direc- 
tions appears only in higher orders in ^(T/fo)". Accord- 
ingly, the right-hand side of Eq. ( [l3| involves only three 
bosons moving in the same direction. This kind of scat- 
tering processes has a divergent rate when the spectrum 
is taken to be strictly linear from the outset, see Eq. ([T]) 
and Fig. [T|a). While Eq. (13) also corresponds to the 
limit of vanishing spectrum nonlinearity, it is crucial that 
the spectrum is linearized after the scattering amplitudes 
are evaluated and the divergences are regularized. 



After integration over 172 and p, Eq. ( 13 ) assumes the 



form 



dt 



dyg{x,y)f{y,t), (14) 



where fix,t) = fqit)\q=2^(^T/hs)x- kernel Q{x,y) is 

given by 



Q{x,y) 



xy{x + y) 



xy{x - y) 



smh\T:{x + y)\ sinh[7r(a; — y)] 
^Ix'ix^ + l)5{x~y), (15) 

D 



and the typical scattering rate is 

Tq-i = 2t: \^K'^s{T/hsf. 



(16) 



The integro-differential equation (14)-(15l can be solved 
exactly. The solution reads [H] 



fix,t)^aoMx)+ dva,ip,{x)e-'^"''^°, (17) 
Jo 
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where rji, = v^iiP' + l)/6 and 

X 



(y5o(x) = v67r — 



sinh(7ra;) ' 
1 



(18) 



3a; 3x 



(2i/2 - l)(5(x- i/) 



sinh[7r(x + z^)] sinh[7r(a; — z^)] 



(19) 



(The singularity in the right-hand side of Eq. ( 19 ) is to 
be understood as the principal value.) The coefficients 
OfQ and OLy in Eq. ( |17p are determined by the initial con- 
ditions, (x^J^ = J^dxip^{x)f{x,0) for /i = 0, i^. 



The first term in the right-hand side of Eq. (171 rep- 



resents a stationary (independent of t) contribution to 
f{x,t). At t oo Eqs. ^ and ^ yield 



hsq/2-lTT' 



(20) 



This result has a clear physical meaning. In general, a 
stationary (equilibrium) solution of the Boltzmann equa- 
tion is not unique. All such solutions 
have the form of the Bose function n 
temperature T. A change of T by 5T generates a cor- 



however, 
q, parametrized by 



rection to N, 



i\t- 



which, to linear order in 5T, indeed 



has the form (|20| with = ^JttJQ{5T /T). On the other 
hand, the energy of the system at i — > cx) coincides with 
that in the initial non-equilibrium state. Thus, the tem- 
perature T characterizing the equilibrium distribution at 
t — )■ cx) is uniquely determined by the initial conditions. 
Choosing Uq as the Bose distribution with this tempera- 
ture, one ensures that ao = in Eq. (17). 



The remaining (time-dependent) term in the right- 
hand side of Eq. (17) describes approach to equilibrium. 
At short times, t <C tq, only the relaxation modes with 
^ {'''o/tY^^ ^ 1 are affected. Since <fu{x) ~ Slx — v) at 



where F{x,t) = f{x,0)e ^ ^/^'^o corresponds to expo- 
nential relaxation with the rate 



127r 



(22) 



The role of the second term in Eq. (21 ) is to ensure the 



energy conservation. The corresponding correction to the 
distribution function [see Eq. ([9])] can be cast in the form 



5Nr, 



dUq 

dT 



ST{t), 6T(t) 



3hs 



dqfq{0)e 



with 1/Tq given by Eq. (22 1. For generic /q(0), the correc- 
tion to temperature ST{t) exhibits non-exponential de- 
pendence on time. 

To summarize, elementary excitations of one- 
dimensional interacting systems are often described in 
the framework of the effective Luttinger liquid theory. 
Both the conventional Luttinger liquid theory '1'3'4| and 
its recent extensions |6l7ll5j provide a set of efficient tools 
for evaluation of various correlation functions. However, 
none of these approaches is capable of describing the ther- 
malization of bosonic quasiparticles because interaction 
between bosons with linear spectrum results in a diver- 
gent inelastic scattering rate. 

In this Letter we demonstrated that the divergences 
are regularized when the nonlinearity of boson spectrum 
is taken into account. We derived and solved the Boltz- 
mann equation describing the fastest equilibration pro- 
cess in the system, namely, thermalization of bosons mov- 
ing in the same direction. The equation describes bosons 
with a linearized (as opposed to strictly linear) spectrum 
and results in a finite relaxation rate that scales with 



1/ » 1, Eq. ([T7| gives f{x,t) (xe ')-*/-ro^ which describes temperature as T^, see Eq. (fiel). Our results are applica 



exponential relaxation with the rate given by the appro- 
priate limit of Eq. ([s]) [g < {T/hs^Y^^, see Eq. 

At t ^ To the high-energy bosons have already re- 
laxed, and thermal bosons (with a; ~ 1 or g ^ T/hs) 
have equilibrated among themselves, although at tem- 
perature that has not yet reached its equilibrium value. 
Indeed, at large t the main contribution to the inte- 
gral in Eq. (17) comes from small v. Approximating 
(p^{x) ~ — 5{x — u) + -^/G/tt (pa{x) and rj^, « v^/& 
find = -f{v,Q), and Eq. (17) yields 



we 



f{x,t)^Fix,t)- ^6/^ipoix) dyF{v,t), (21) 



ble to both fermions and bosons with strong long-range 
repulsion. 
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In Sec. [T] of this Supplemental Material we derive the 
dispersion relation for phonons in the Wigner crystal. 
The derivation of the on-shell amplitude of the lead- 
ing inelastic scattering process is discussed in detail in 
Sec. [2] In particular, in Sec. |2.1| we address the screened 
Coulomb interaction, and in Sec. 2.2 we demonstrate 
that the amplitude vanishes for integrable models with 
V{x) oc l/sinh^(cpa:) and V{x) oc l/x^. Finally, in Sec.js] 
we present the exact solution of the linearized Boltzmann 
equation. 



1. PHONONS IN A WIGNER CRYSTAL 

We consider a system of N identical spinless particles 
of mass m described by the Hamiltonian [sec Eq. (4) of 
the Letter! 



2m 2 



(1.1) 



where pi and xi are, respectively, the momentum and 
position of the Zth particle {I — \, . . . , N), and V{x) is the 
interaction potential. Expanding the potential energy in 
Eq. (1.1) to leading order in \ui — we obtain the 



Hamiltonian of a harmonic chain, 



2m ' 4 ■ 
I 1,1' 



where we introduced the notation 
_ d^Vjx) 



dx^' 



(1.2) 



(1.3) 



Sqq'. Substitution into Eq. (1.2) then yields 



where the phonon frequencies ujq are given by 



(1.6) 



(1.7) 



1=1 



Since Eq. (1.71 does not contain h, the same result can 



be obtained by solving classical equations of motion. 

At small q and for V{x) decaying with the distance 
faster than 1/x^, Eq. (1.7) reduces to Eq. (3) of the Let- 
ter, 



ujq = s\q\(l - ^q'^ 



with 



where 



Vrr. 



24 V22 



(1.8) 



(1.9) 



(1.10) 



1=1 



x=l I p 



Modifications of Eq. (1.8) for potentials decaying with 
the distance as Xjx^ (screened Coulomb potential) and 
Xjx^ (Calogero-Sutherland model) are discussed below in 
Sec. |2.1|and Sec. |2.2[ respectively. 



It is convenient to write ui and p; in the second-quantized 
representation. 



SCATTERING AMPLITUDE 



Pi 



E 



2mNujq ' 



{bq + blq)e^''\ (1.4) 



where the phonon creation and annihilation operators 
satisfy the canonical commutation relation — 



We start with the general expression for the on-shell 
scattering amplitude in order H"^ derived in Ref. [T], 



A 



(2.1) 



where 



2 



A 



h{<ii,q2)h{q'i,q'2) , h{<i2,-q'i)h{qi,-q2) , h{qi,-q'i)h{q2,-q'2) , -m 



2 A2 + y/4(9l,'Z2,-«;) 



with 



C30 

h{qi-,q2) = ^y/^^|sin[(qi+g2)^]-sin(giO-sin((j20}, 

00 

/4(gi, 92, 93) = X! ^/^^^{^ ~ cos(giO - cos(g20 - cosi^q^l) - cos[(qi + 92 + 93)^] 

+ cos[(gi + 52)'] + cos[(g2 + q:^)l] + cos[(qi + q3)Z]|. 



(2.2) 

(2.3) 
(2.4) 



Expanding Eqs. (2.2|-(2.4| to first order in qi,q2,q'i,q'2, and taking into account the momentum and energy conser- 
vation laws, qi + q2 = q'l + 92 ^■'^d + ojq^ — ujq'^ + ujq'^, wc arrivc at Eq. (5) of the Letter, 



with 



Here 



A = 



V22V44 — V33 



1/2 



919291 92 



1/2 



Ap^Vi^ l&pWl g^g'^O 1 F[Aq] 



^'g+g' - ^g " ^g' 4py22 



^33 -^[A,] J ■ 



^ ^ dq mp 



^y/^)/[l-cos(gO], 



(2.5) 



(2.6) 



(2.7) 



and the functional J- is defined as 

/(9 + 9') f{q) 



Hf{q)] 



9' 



/(9O 
9' ■ 



(2.8) 



Interestingly, although the right-hand side of Eq. (2.6) 



depends on the functional form of the interaction poten- 
tial, it is independent of the interaction strength: mul- 
tiplication of V{x) by an arbitrary constant leaves A in- 
variant. 

Equation (2.61 is valid provided that l^m a-i'e finite, 
i.e., for interaction potentials decaying with the dis- 
tance faster than 1/x. Further simplification is pos- 
sible if Vn,n+2 [see Eq. (1.10)] are finite as well, i.e., 
for lim^^^oo x^V{x) — 0. Expanding A{q) and -6(9) in 



Eq. (2.6) to fourth order in g, we find 



A, = sV(l-2e9') 



mp 



33 



(2.9) 



into Eq. (2.6) results in 



with s and ^ given by Eq. (1.9). Substitution of Eq. (2.9) 



A = 



V24V44 — V33V35 

4^2^221^24 



(2.10) 



2.1. Screened Coulomb interaction 

We turn now to the important in practice case of the 
Coulomb potential screened by a remote gate at distance 
d from the Wigner crystal. 



V{x) = 



Va;^ + 4^2 



(2.11) 



For this potential Vnn are finite and determined by the 
behavior of the potential at short distances x < d. For 
pd^ 1 we find 



Vr,„ = {-irnleY+Hn{pd). 



(2.12) 



Vn,n+2, however, diverge due to slow decay of V{x) at 

Vn.n+2 = (-l)"(n + 2)! e2d2^n+35^ (2.13) 

where S = ^ logarithmically divergent sum. 



(In writing Eq. (2.13), we neglected the contribution from 



I < pd.) Treating 5 as a large, but finite number (or, 
equivalently, introducing a cutoff I < Iq with ^ pd)^ 
we find V35/V24 — —bp. Substituting this relation into 
Eq. (2.10) and taking into account Eq. (2.12), we obtain 



A 



-3/4. 



(2.14) 
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An accurate derivation of Eq. (2.14) should rely on for n = and n 



Eq. (2.6) rather than Eq. (2.10). Even though V24 di 



verges, the dispersion relation can still be cast in the 
form of Eq. (1.7), but with g-dependent ^. For small q, 



q < 1/pd < 1, 



we find 



59(1 - Lq 



with 



v/(2eV/TO) Hpd) , 



2 \n{pd) 



(2.15) 



ln(l/g). (2.16) 



In the same limit, the functions Aq and i?^, see Eq. (2.7 1, 
are given by 

Aq^sY{l~2£,qq^), = -sV(3-5e,g'). (2.17) 



Substituting Eq. (2.17) into Eq. (2.6 1 and taking into 
account Eq. (2.12), we recover Eq. (2.14). 



2.2. Integrable models 

In integrable models |2| conservation laws forbid a re- 
distribution of momenta between colliding particles. In 
other words, collisions do not lead to relaxation in such 
models. This lack of relaxation should not depend on 
whether the system is described in terms of the bare par- 
ticles or in terms of the quasiparticles. In accordance 
with this observation, the amplitude (2.1 ), as well as the 



amplitudes of higher-order scattering processes, are ex- 
pected to vanish for all integrable models. This vanishing 
of the scattering amplitude also serves as an independent 
check of the validity of Eqs. (2.6 1 and (2.10) above. 
Consider the potential [5] 



sinh {cpx) 



(2.18) 



For this potential 



— V{x)^d'' — V{x)x- 



and Eqs. ( 1.3 1 and ( 1.10 ) yield 



Another useful relation, 



1 d 

^°^-"2d^^°°' 
can be derived by comparing the expansions 



(2.19) 



(2.20) 



00 00 



-Iclrfi 



( |2.20[ ), we find 



2. Using the identities (2.19) and 



= -2 



{cpy 



5c" ' 



which gives 
V24V44 



{cpf d^V^a d^Vo^ 
2 dc^ dd^ ' 



(2.21) 



(2.22) 



Substitution into Eq. (2.10) then yields A = 0. 

At c ^ 1 and in the Wigner crystal limit Vq ^ 



T?(?e^'^ jmi? ^ the potential (2.181 realizes the Toda lat- 
tice model. In this regime 



V{x) = 4Vbe-2'=''l^l 



(2.23) 



and all but / = 1 contributions to the sum in Eq. (1.10) 



can be neglected (this corresponds to the nearest neigh- 
bors interaction in the Wigner crystal). With these ap- 
proximations. 



K„„-4yo(-l)"(2cp)"e-2-, 
and the relation V24V44 = V33V35 is obvious. 



At c ^ 1 the potential (2.18) becomes 



V{x) 



(2.24) 



(2.25) 



with a = Vo{cp)~^ {a ^ It? /m in the Wigner crystal 
regime), which corresponds to the integrable Calogero- 
Sutherland model [2]. For this potential Vnn are finite, 



Vnn^{-lT{n + l)\{^^/Q)ap 



n+2 



(2.26) 



but Vn n+2 diverge, and Eq. (2.10) is inapplicable 



Equation ( 2.18 1 can be viewed as a version of Eq. ( 2.25 ) 
with a long-distance cutoff. Importantly, for the 
Calogero-Sutherland model the sum entering Vn.n+2 [see 
Eq. ( |1.10 )] diverges as a power-law rather than logarith- 
mically, as is the case for the screened Coulomb interac- 
tion, see Sec. 2T Accordingly, the value of A depends 
on the cutoff scheme and is not universal. It is therefore 
important to demonstrate vanishing of the scattering am- 
plitude directly for the Calogero-Sutherland model (2.25) 
instead of relying on c ^ limit of Eq. (2.181. 



For the potential (2.25) we find 



A„ 



1 



Br, 



2tt 



2 4 



(2.27) 



Substituting Eq. (2.27) into Eq. (2.6) and using 



Eq. (2.26), we indeed find A = 0. Note that for the 



Calogero-Sutherland model 



= ^q^sq{l-q/2Tlf 



(2.28) 



i.e., the nonlinear term in the dispersion relation is 
quadratic in q rather than cubic as in Eq. (1.8 1; this 



; = 1 r?j.= l 



property is specific [2] for the inverse-square interaction 
potential. 
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3. SOLUTION OF THE LINEARIZED 
BOLTZMANN EQUATION 

Solutions of the equation [see Eqs. (14)-(15) of the 
Letter] 



d_ 
di 



Q{x,y) = 



'0 / dyGix,y)fiy,t). 



(3.1) 



xy{x + y) 



xy{x - y) 



sinh[7r(a; + y)] sinh[7r(a; — y)] 

+ lx\x^ + l)S{x-y) (3.2) 
6 

have the form f{x,t) = </?(a;) exp(— Tyt/rg), where ip{x) 
satisfies the integral equation 



r]ip{x) = / dyg{x,y)ip{y). 
Jo 



(3.3) 



Remarkably, the eigenvalue problem ( |3.3[ ) can be 
solved exactly. First, we formally extend f{x) to neg- 
ative x according to 

ipi-x) = ip{x) (3.4) 

and rewrite Eq. (3.3 1 as 

^ ^ 1 2. 2 , 1^ . N , f°°, xy{x + y) 

6 smh|_7r(.T + y)J 

(3.5) 



Next, we multiply both sides of Eq. (3.5) by e*'-^ and 



integrate over x. This transforms the integral equation 
(3.5) into a differential equation, which can be written as 



(/i2 + /i^)^(C) = 6r;^(C). 
Here <^(C) is the Fourier transform of ip{x), 

^(C)-^(-C)= /dxe^Xx), 



and the operator ft,^ is given by 

3 

k = - 



2cosh^(C/2) 



(3.6) 



(3.7) 



(3.8) 



Equation ( 3.8 1 coincides with the Hamiltonian of a par- 



ticle moving in one dimension in the presence of the re- 
flectionless Poschl- Teller potential. The corresponding 
eigenvalue problem, 



is discussed in detail in, e.g., [3]. Even-parity eigenstates 
of h(^, which are of interest here [see Eq. (3.7)], include 
the ground state with the eigenvalue e = — 1, 



V'-i(C) 



1 



cosh^(C/2) 



(3.10) 



and a continuum of states with eigenvalues e — v^, 

MO = (- I + t^-h 2)(- I + 2 '^"^ 2) '"'("^^ 

(3.11) 

Obviously, the eigenstates of h(^, s ee Eq. ( |3.9[ ), are also 



eigenstates of hj^ + see Eq. (|3.6|). The corresponding 



eigenvalues are related according to 77 = e(e+l)/6, which 
gives 770 = for the bound state (e = ^1), and r/^ = 
v'^{v'^ + l)/6 for the continuum (e — v^). 

Carrying out the inverse Fourier transform of Eqs. 
([3I0I) and ([3ll|), we find Eqs. (18) and (19) of the Letter, 



ipo[x) = V67r 
ipt,{x) 



sinh(7ra;) 
1 



(3.12) 



V(!y2 + l)(4i.2 + i) 
3a; 



(2j/2 - l)5{x-v) 



3a; 



+ . , , , rr + . , , , rT)>, (3.13) 

sinh[7r(a; + v)\ sinh[7r(a; — v)\ 



with the singularity in the right-hand side of Eq. (3.13) 
understood as the principal value. 



The eigenfunctions (3.12) and (3.13) are normalized 
according to 



dxLpQ{x) — l, I dxipu{x)(p^'{x) — 5{v ~ v') 
'0 Jq 

(3.14) 

[ipo{x) is orthogonal to ipuix) for any v], and form a com- 
plete set, 



(poix)tpo{y) + dvtp„{x)tpu{y) = 5{x- y). (3.15) 

^0 



Thus, the general solution of Eqs. (3.1)-(3.2) can be 
written as an expansion in tpo and ipi,, 

f{x,t)^ao'Poix)+ / d7^a^V^(a;)e-''"*/^«, (3.16) 



kMC) = eMO, 



(3.9) see Eq. (17) of the Letter. 
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